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shells. The balance between self and exterior doping represents the best strategy for 
boron buckyball synthesis. 
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I. INTRODUCTION 



Boron analogues^ to carbon sp^ structures usually take the form of compounds where 
an other element acts as an electron donor. In these structures, the B~ anions behave 
like carbon atoms. This is the case of MgB2'^, the boron analogue of graphite. For less 
extended structures this purely exterior doping can be unfeasible and impractical, e.g. giving 
60 electrons to a boron buckyball is no simple task. Fortunately, simulations have shown 
that pure boron systems can partially overcome the low electron count on their own through 
self-doping. This mechanism explains the stability of small planar sp^ clusterd^'^, the 
relative stability of the all boron a-sheet^ and other structure^^^^. In contrast with carbon, 
these all-boron structures were also expected to display polymorphismpHH ^^^i^ multi-center 
bondin^SEl allowing new behaviors for well-known structures. On this basis, the Bgo 
fuUerene has attracted a lot of attentioiP^^. 

However, no boron fuUerene were observed experimentally. Consequently, it was shown 
that other less symmetric bulk-like precursors, dubbed core-shell structures, are energetically 
preferred ^^ * ^^^ . Self-doping is thus insufficient to guarantee sp^ bonding in extended boron. 
Nevertheless, it can be argued that synthesis of these systems remains possible if both self 
and exterior doping are present. Endohedral metalloborofuUereneJ^SESl were found to be 
energetically more stable than the corresponding core-shell cluster^, suggesting a possible 
synthesis pathway which is very similar to the one used in the synthesis of endohedral carbon 
fuUerenes. The recent observation of boron nanotubes inside Mg rich catalytic pores^'' or on 
catalytic surfaced give even more weight to this scenario. 

Yet, some problems would remain. The polymorphism of boron can also prove to be a 
problem since there would be no energetically motivated driving force towards some well- 
defined structur^. A successful synthesis would lead to disordered structures which could 
negatively impact their characterization and properties. Furthermore, polymorphism is also 
related to high reactivity: endohedral boron fuUerenes could still have the possibility of 
forming outward sp^ boron structures by merging with other boron clusters during synthesis. 

In our paper, we offer a critical reassessment of the role of doping in the boron fuUerenes 
and show how it can overcome the above mentioned problems. We first rationalize the role 
of self-doping in the pure Bgo borofuUerene. The different Bgo fuUerene isomers are formed 
of a Bgo backbone with an additional 20 filling boron atoms, which can be placed in the 
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FIG. 1. Wannier functions representing the bonding pattern of a filled pentagon for the C3 17-3 
isomer. The 3c-2e Wannier functions are represented in red and the 4c-2e Wannier functions are 
in blue. The neighboring filled hexagons are shaded to clarify the local environment. A schematic 
representation of this pentagonal environment is shown at the bottom right. This coloring scheme 
is reused in Figure |2j 



center of hexagons or pentagons. These boron atoms act as dopants. We show that filled 
pentagons contribute to lower the global strain of the fuUerene, but a specific type of local 
environment is necessary to satisfy the electron-counting issues arising from the self-doping 
rule for filled pentagons. Additionally, it will be shown that delocalized tt bond percolation 
(partial aromaticity) plays a stabilizing role in these isomers. All of these findings support 
the fact that, for these open-shell clusters, the pentagonal pyramids represent the crucial 
building block for the stability of the Bso fuUerenes. This leads us to formulate a generalized 
isolated filled pentagon rule (IPFR). Having understood the consequences of the strain and 
of self-doping in pure borofuUerenes, we then focus on the role of exterior doping, i.e. charge 
transfer via endohedral centers (seeds) or anions. We show that exterior doping not only 
further stabilizes these clusters but also diminishes their reactivities. Furthermore, the 
overall polymorphism can be broken by the interaction with a seed. Hence, endohedral 
metalloborofuUerenes possess all the required properties to overcome the known hurdles 
inhibiting borofuUerene synthesis. 
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FIG. 2. The different B^o isomers with increasing number of fihed pentagons. They were chosen 
to exemphfy the role of the filled pentagonal unit inside these structures. The different local 
environments of the filled pentagons which are present in each isomer are sketched below them 
along with their frequency. The colored sections of these sketches represent a neighboring filled 
hexagon. Also displayed is the energy difference ( AE in meV) between the isomer and the Ih 20-0, 
its HOMO-LUMO gap (Eg in eV) and its average pyramidalization angle (6). 



II. METHODS 



We use the Kohn-Sham approach to density-functional theory (DFT) as implemented in 
the BigDFT code^, which uses a systematic real-space wavelet basis. We considered the 
basis converged when an accuracy of 0.5 meV/atom was reached for the total energy and 
1 meV/A for the forces. This corresponds to a uniform wavelet grid with a spacing of 0.4 
bohrs with an extent of 11.5 bohrs for the coarse grid and 2.9 bohrs for the fine grid. The 
structural relaxations used the FIRE algorithm^'' and were stopped when the forces were 
below 5 meV/A. In addition HGH pseudopotentials^^ in the Krack variant^ were used with 
isolated boundary conditions. Hence, no supercell approximation was needed. 

Our calculations used the GGA-PBE^^ exchange-correlation functional, which was shown 
to give the same energy organization as coupled cluster methods for small boron clusters 
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(in contrast to other hybrid functionals)^. The GGA-PBE was also shown to yield more 
accurate results than LDA in these systems, even if some deviation from diffusion quantum 
Monte Carlo was seeii^. 

The Wannier functions were extracted from the Kohn-Sham wavefunctions using the 
Wannier90 cod^^. They were considered converged when the total spread varied by less 
then 10~^^A^. An initial projection on sp^ functions located on the atoms and with the 
z axis pointing along the radial direction was needed. In order to quantify the degree of 
strain inherent to each isomer, we have used the Pi Orbital Axis Vector (P0AV2) method^ 
developed for carbon fuUerenes. The P0AV2 analysis of the structure was also done with 
a code developed for this purpose. The pyramidalization angle is defined as the mean 
deviation from perfect orthogonality between the a and tt orbitals. The multi-center bonding 
revealed by the Wannier analysis (or other local orbitals analysid^ coupled to the high 
coordination number of boron atoms (5 or 6) adds some difficulties in the POAV analysis 
of these systems. Luckily, a study of the Wannier functions revealed that the filling atoms 
do not possess local orbitals with tt^ character. They thus literally behave as doping atoms 
which contributes a electrons in three 3c-2c bonds. Hence, there is no point in evaluating 
there tt misalignment. On the other hand, the backbone atoms do possess tt^ character and 
the P0AV2 analysis was done only for these atoms. The still relatively high coordination 
of these atoms (5) did not pose a problem because the 3c-2e Wannier centers were located 
on the plane defined by the corresponding three atoms and thus using any of these atoms 
yields the same pyramidalization angle. Thus, for simplicity, only the backbone atoms were 
used for this analysis which reduced the effective coordination to three. 

In order to be concise, we do not discuss the pseudo-Jahn- Teller distortiori^ which was 
found to yield a slightly more stable T^ 20-0 fuUerene. This distortion has no impact on 
the conclusion of this work. The associated puckering of the filling atoms does not impact 
the global pyramidalization angle because they are alternatively moved inside and outside 
the structure. Alternatively, this distortion can be understood in the IFPR since for the 
20-0 structure the filling atoms are the most reactive sites. In the rest of the text, the 
fuUerenes are referenced according to their symmetry and their number of filled hexagons 
and pentagon^. For example, the first proposed fuUeren^^ becomes the Ih 20-0 and the 
volleyball structure^^ becomes the T^ 8-12. 
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III. RESULTS AND DISCUSSION 



A. Effect of filled pentagons on pure borofuUerenes 

In the fuUerene structure, the role of the pentagons is to aUow for curvature by increasing 
the pyramidahzation of its cornerd^S. Hence, most of the strain is stored in the pentagons. 
In carbon, this leads to the isolated pentagon rule (IPRpOlMJ which is a rather strong re- 
quirement on the stability of pure neutral carbon fuUerenes. In the boron fuUerenes, the 
pentagons essentially play the same role. The only caveat is that the filling boron atoms 
can be placed to lower the global strain of the fuUerene by increasing their pyramidahzation 
angle. Since these atoms serve mostly as dopants they are not constrained to a sp^ bonding 
pattern, in a way reminiscent of the carbon atoms in 77^-complexes^^ of carbon fuUerenes, 
and they do not contribute to the strain. The strain is characterized by the average pyrami- 
dahzation angle (see Methods) and is presented in Figure [2| As can be seen, the inclusion of 
filled pentagons always releases strain and the energetically preferred structure (Dad 14-6) 
sports the smallest angle of all the isomers. 

Nevertheless, not all the isomers are stabilized when compared to the Ih 20-0. Indeed, 
the position of the 20 filling atoms with respect to the 60 atoms backbone cage appears to 
be crucial. They should be distributed such that the self-doping rule for boron sp^ systems 
is satisfied: filled rings should be bordered by alternating filled and empty neighbor rings. 
Of course, this cannot be realized for the filled pentagons and they will thus inevitably 
be frustrated. The most stable local environment for filled pentagons, presented in Figure 
[T| is the closest one can find to the self-doping requirement: almost all backbone atoms 
contribute half their electrons to the pentagon, except the atom located at the junction of 
the two filled neighbors which contributes only one. This atom has a coordination number 
of six, instead of a coordination of five which is normally found in the backbone. It thus 
behaves more like a dopant atom than like a backbone atom and we can see this side of the 
pentagon has three inter-penetrating fiUed-hexagons (see supplementary information). It is 
this boron atom which sports the largest pyramidahzation angle but because of its almost 
dopant nature, it does not contribute to the strain. 

These two considerations (strain and self-doping) explain the total energy diflFerences 
observed in Figure [2} The inclusion of only one or two filled pentagons always results in a 
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Dsd 14-6 




FIG. 3. Schlegel diagrams with superposed Wannier states for the Ih 20-0 and both Csy 14-6 
structures. The 3c-2e bonds are in red and the 4c-2e bonds are in blue. The delocahzed tTz 
bonds are displayed in green on the right panels. For clarity, the filling atoms are shown as big 
purple circles. The Ih 20-0 is formed of alternating 3c-2e bonds and 4c-2e bonds in the same 
Kekule structure then the one observed in Ceo- In the case of the 14-6, we can see a complete 
percolation of the tt^ states on the equator. For the D^d 14-6 (B), the tTz states span out of the 
poles in three lightning shaped arms that do not meet. 

destabilized fuUerene because the self-doping rule cannot be achieved. The first stabilized 
fuUerene is the C3 17-3, where the three-filled pentagons are clustered around one hexagon 
which becomes the pole of symmetry. All other stabilized structures, are completely formed 
of pentagons with the same local environment, except the 10-10 which also possess other 
types of pentagons but which are placed to complement each other. This demonstrates 
that such an environment is the most stable for the Bgo and that the increased stability 
in these isomers stems mostly by global release of strain as denoted by their decreased 
average pyramidalization angles. This balance between strain and electronic issues form a 
generalization of the IPR, which we call the isolated filled pentagons (IFPR) rule, in the 
borofuUerenes. 
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B. Partial aromaticity and reactivity 



To quantify the sp^ character of the cages, we perform a decomposition of its orbitals 
in terms of locahzed bonds via the Wannier transformatioip21 (^see Methods section). The 
Kekule bonding pattern found in the original Ih 20-0 is similar to the one found in Ceo except 
for the multi-center character of the bonds: there is a double four-center-two-electron (4c-2e) 
cr-bond located between two adjacent fiUed-hexagons and a three-center- two-electron (3c-2e) 
cr-bond located on the edge of the empty pentagons, as seen on the Schlegel diagrams of 
Figure [sj This is also similar to the bonding pattern found in the boron a-sheet^^, with 
one notable exception: the tt^ states are now partly localized inside the double 4c-2e bonds. 
This lowering of the aromaticity^^ further destabilizes fuUerene structures with respect to 
their planar counterparts. 

For the other isomers, the bonding pattern of the filled pentagons slightly changes because 
of the differing symmetry, i.e. not all surrounding hexagons are necessarily equivalent. 
Nevertheless, a Wannier analysis shows that most states remain the same. Notable changes 
stems from local release of strain which are associated to the presence of strongly delocalized 
Tiz 7c-2e states. As can be seen in Figure |3| in the Dad 14-6 these delocalized tt^ electrons 
percolate around the equator of the isomer. This equatorial aromaticity is analogous to the 
one observed in C70, albeit weaker : the equatorial angles are of 10.5° for this isomer while 
they are of S.TS^I^for Ctq. Consequently, the atoms of the equator are thus expected to be 
less reactive than their 20-0 counterparts. At the same time, other sections of this isomer 
seem far more prone to reactions. The chemical hardness, determined by the HOMO-LUMO 
gap, can be used on the DFT level as a good representation of the reactivity trends in these 
isomers. Furthermore, the local strain coupled with the atomic contributions to frontier 
orbitals was shown to yield accurate reactivity maps for convex clusters^^. The local strain 
on the backbone atoms of the poles (11.92°) is larger than for the Ih 20-0 (11.46°) which 
leads to a higher strain activated reactivity. The non-degenerate LUMO level of this isomer 
is also located near these atoms suggesting that they indeed represents the preferred sites 
for reduction reactions. 

One can thus wonder if a different distribution of the same number of filled pentagons 
can lead to a decreased reactivity while maintaining the global release of strain. The IFPR 
suggests that this is not possible because release of global strain is associated with its 
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ScsNOBso 


^80^ 


AE [eV] 


Eg [eV] 


AE [eV] 


Eg [eV] 


Ih 20-0 


-10.859 


0.592 


-3.563 


0.132 


Dad 14-6 


-10.973 


0.229 


-4.591 


0.550 


Dsd 14-6(B) 


-12.976 


0.385 


-4.286 


0.175 



TABLE I. Binding energies and HOMO-LUMO separation for the endohedral Sc3N@B8o fullerene 
and the B^q^ ion. Binding energies are mesured with respect to the energy of the pristine cage (plus 
the energy of the added cluster). 

local increase in some zone of the structure. Thus increased stability comes at the price 
of increased reactivity. This is demonstrated within a second 14-6(B) isomer. This 
structure is formed with exactly the same pentagonal environment but the filled pentagons 
now lie on the equator. This corresponds to the maximally separated uniform distribution 
of filled pentagons. In this case, both the global pyramidalization angle (11.35°) and the 
minimal angle (10.60°) are higher than the ones of the 14-6. Now, its HOMO-LUMO 
gap is almost five times higher suggesting a much lower reactivity. An inspection of its 
doubly degenerated LUMO reveals that it is now delocalized on the equator, suggesting 
that it possesses less specific reactive sites. The same kind of delocalized tt^ electrons is 
seen in this isomer, but because of the new pentagonal arrangement, they are this time 
arching from the poles in three lightning shaped arms that terminate on the filled pentagons 
(represented in the lower right section of Figure [3]) . This points to a lower reactivity of 
the poles. These tt^ states do not percolate leading to a lower resonance energy then the 
14-6. The increased strain and reduced resonance energy explain the energy diflFerence 
between the two 14-6 isomers. 

C. Exterior v.s. self-doping 

Having identified the consequences of strain and self-doping, we can now address the 
infiuence of exterior doping, i.e. the interaction with non boron atoms. As mentioned earlier, 
metallic clusters are envisioned to act as growth centers for endohedral metalloborofuUerenes. 
The induced stabilization of each endohedral isomer will depend greatly on the nature, 
symmetry and orientation of the seed. Thus, the identification of a good candidate for 
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synthesis would necessitate a complete zoology of metallic clusters coupled with their effect 
on different isomers. This is beyond the scope of this paper. We will rather focus on the 
general effects of the inclusion of the metallic clusters on the pentagonal unit, by considering 
two basic examples: the SC3N cluster and the B^q anions. 

As can be seen in Table [T| all isomers are stabilized by the inclusion of the SC3N "seed" . 

For the Ih 20-0, we find that the nitrogen atom departs from the plane formed by the Sc 
to form a covalent bond with a filling boron atom, thus reproducing the results of Peng et 
al.'^. This is consistent with the reactivity picture described by the LUMO and the strain. A 
Bader analysis of the electron density reveals that the Sc atoms have lost 1.5 electrons while 
the N atom gained 2.5 electrons. This corresponds to a net transfer of 2 electrons to the 
boron system. A Jahn- Teller deformation is observed at this filling since the LUMO of the 
Ih 20-0 is threefold degenerate. This reduces the symmetry of the system to Ci leading to a 
decrease of the HOMO-LUMO gap of the system to 0.592 eV. The global pyramidalization 
angle remained the same (11.46°) showing that the endohedral cluster in the Ih 20-0 do not 
impact the global strain. Since the symmetry is decreased, there is now a greater scatter in 
the local strain leading to some preferred reactive sites which are located close to the N-B 
bond. 

In the Dad 14-6 isomers, things are quite different. The metallic cluster adopts a quasi- 
planar structure which places itself on the plane containing the equator. This shows that 
no covalent bonds between the cluster and the fuUerenes are formed, which is confirmed by 
the electronic density. This also agrees with the previous reactivity picture: the LUMO of 
these systems is formed of backbone tt^ states. Interior covalent bonding with these atoms 
is therefore unlikely: it inevitably leads to an increase of the global strain which greatly 
destabilizes these states. This is the origin of the "topological" protection invoked in the 
literatur^^. Nevertheless, the global pyramidalization angle of the 14-6(B) was greatly 
reduced (11.09°) and is now similar with the Dad 14-6. Consequently, as we can see in Table 
[l| the stability of the Dad 14-6(B) is greatly increased. It becomes the most energetically 
stable isomer. This can be traced to electrostatic effects. A Bader analysis of the charge 
density displays a 3.2 electron transfer from the cluster to the cages (with +1.6e Sc atoms 
and -1.6e N atom). These extra electrons are mostly clustered around the equator which 
maximizes the electrostatic interaction with the positively charged Sc atoms. Its increased 
stability is thus a consequence of the planar geometry of the seed. By maximizing the 
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electrostatic interactions, by varying the shape and the symmetry of the metaUic clusters, 
it is thus possible to maximize the driving force which will select a given fuUerene structure 
during growth. 

Interestingly, the seed increases the HOMO-LUMO gap of the D^d 14-6 suggesting a 
decrease in reactivity. This can be maximized if the metallic cluster is chosen to yield an 
electron count that matches the distribution of the electronic levels of the neutral fuUerene. 
To this end, we have studied the corresponding B^^ ions. Here, we see that the stabilization 
caused by the population of the tt^ electrons is maximal for the D^d 14-6 which stays the 
most energetically favored isomer. Since the LUMO of this isomer is not degenerate, there 
is no Jahn- Teller deformation in this system and the pyramidalization angles stays the 
same. The rigid band approximation works well in these systems, since a difference of the 
anion and neutral densities show that the extra charge is located on the LUMOs and the 
energy levels are basically unperturbed. The extra electrons are thus distributed in the tt 
states of the backbone atoms of the poles of the D^d 14-6 and it is the increased resonance 
energy which is responsible for the stabilization and decreased reactivity of the structure. 
Unsurprisingly, in the case of the D^d 14-6(B), there is a Jahn- Teller deformation since in 
its neutral state the LUMO is twice degenerated. This deformation allows it to decrease its 
global pyramidalization angle (11.27°) which is responsible for its further stabilization. 

IV. SUMMARY AND CONCLUSIONS 

To rationalize the effect of exterior and self-doping, we carried out simulations on the 
Bso isomers. This has lead us to formulate the generalized isolated filled pentagon rule 
(IFPR) which states that a peculiar environment is necessary for the stability of the filled 
pentagons. Because of polymorphism, the IFPR is not a strong constraint on stability, as 
the IPR, but rather a guideline for constructing the most stable known structures. It is also 
shown that increased tt conjugation enabled by local decrease of strain plays a large role in 
the stabilization of certain isomers. The IFPR isomers are thus reminiscent of C70 where 
equatorial aromaticity plays a stronger role. This is in good agreement with the simulated 
ring currents for the Ih 20-0^. Further study of the effect of filled pentagons on the ring 
currents of Bgo would be needed. Furthermore, the IFPR clearly indicates that the synthesis 
of pure borofuUerenes faces difficult challenges because increased stability comes at the price 
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of increased strain induced reactivity. 

The IFPR should be vahd for larger boron fuUerenes. In this case, the increase of the 
hexagon to pentagon ratio enables greater freedom to determine the local environment of 
the filled pentagons. It then becomes possible to maintain a same filling pattern than the 
a-sheet, or other stable sheets, for the hexagons located far from the pentagons all the 
while modifying the ones close to the pentagons to maximize stability. However, the relative 
energy differences should drop as the size of the fuUerene increases since the effect of the 
curvature diminishes. Larger fuUerenes are thus expected to have a denser energy spectrum 
recovering the stronger polymorphisnJ^ of the boron sheets. 

Self-doping alone is not enough to guarantee the existence of the fuUerene structures. It is 
thus necessary to explore other means of stabilizing these structures. The use of endohedral 
"seeds" for cluster growth could provide further doping leading to greater stabilization. In 
an effort to understand if a future synthesis of endohedral borofuUerenes is feasible, we have 
studied the effect of this exterior doping on the IFPR. Of course, we have restricted ourself 
to IPR fuUerenes and one must keep in mind that in endohedral carbon fuUerenes departure 
from the IPR is well-knownP' and that this can also be the case in endohedral borofuUerenes. 
In the two cases considered, namely the SC3N cluster and the B^^ anions, the stability of the 
isomers were systematically increased. The effect of the IFPR was found to inhibit direct 
reaction with the endohedral cluster since they would inevitably increase the strain of the 
fuUerene shell. Furthermore, it is shown, that during growth, the symmetry of the metallic 
cluster will select a given isomer which will minimize the electrostatic effects. In other words, 
the "seed" can restore a greater energetically motivated ( 2 eV) driving force. We have also 
demonstrated that electron addition can decrease the reactivity of the borofuUerenes thus 
reducing the probability of forming outwards sp^ structures during synthesis. We thus show 
that a metallic seed has all the requirements needed to solve the difficulties restricting the 
synthesis of boron fuUerenes. 
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